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1 Introduction 



In the frame of the Standard Model the mass eigenstates of the quark fields are not 
eigenstates of the weak interactions: the two set of bases are connected by a unitary 
transformation which is represented by a three times three unitary matrix. This fact 
accounts for flavour changing charged currents and flavour changing neutral currents The 
mechanism of flavour mixing was originally proposed to account for the different ampli- 
tudes of the decays of the muon and of the down and strange quarks, all mediated by 
weak charged currents Jl||. It was extended to three quark generations in 1973 p| (the 
beauty quark, down type member of the third generation was later discovered in 1977) as 
a possible description of the CP violation in K° decays within the frame of the Standard 
Model without introducing new (so called Super Weak) interactions. Once removed all 
arbitrary phases, the three by three unitary matrix must contain four independent pa- 
rameters, three rotation angles in the quark field space and a phase, which introduces an 
imaginary part in the Hamiltonian. Therefore the amplitudes for a process and for its CP 
conjugated (obtained by hermitian conjugation) may differ. 

Several representations of the the Cabibbo-Kobayashi-Maskawa Jl|, |2|] mixing matrix 
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were proposed in the past. The Wolfenstein parametrisation || is adopted in this paper as 
it naturally describes the measured hierarchy among the parameters. The matrix elements 
are expressed in terms of the four parameters A, A, p and r\ 

( 1 - A 2 /2 A AX 3 {p - irj) \ 

V C km= -A l-A 2 /2 AX 2 \+0(X 4 ). (2) 

V AA 3 (1 -p-irj) -AX 2 1 / 

A simple extension valid to (9(A 6 ) is described e.g. in (||. The parameter r\ is the 
complex phase accounting for CP violation of the weak interactions. The A and A pa- 
rameters are known with a good accuracy (~ 1% and ~ 4%, respectively), while many 
contributions to extract p and r\ from the available measurements exist in the literature 
H, H |[. To this purpose, the measurements of the CP violation parameter in neutral 
Kaon decay je^l, of the difference between the mass eigenstates in the B\ — B^ system 

and the lower limit on the difference between the mass eigen- 



Arrid and of the ratio , 

states in the B° s — B® system Am s can be used. On the other hand, since the only direct 
experimental evidence for CP violation is given by the fact that \e^\ ^ and the effect 
could be explained in term of models proposed in alternative to the Standard Model (see, 
for instance, J?], U and references therein), it is suggested to remove the constraint coming 
from the neutral Kaon system and to investigate the results on parameter rj or to test the 
hypothesis of a real Vckm matrix. In || two different procedures have been exploited. In 
the first one the two parameters p and 7] has been fitted to the experimental values of all 
available constraints described therein, apart that coming from measurement of neutral 
kaon mixing. In the second one the same fit has been computed but forcing rj to zero. 
The results of the two procedures in || are opposite. In this letter the second procedure 
has been followed but with a different statistical approach, the Best Linear Unbiased Es- 
timator P], and the conclusion is opposite to the second procedure in H but in agreement 
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with the first one. In this letter, the inclusion of different data sets is also discussed and 
the corresponding results are presented. 



2 Measurements and constraints on Vckm parame- 
ters 

The A parameter is the sine of the Cabibbo angle |L(| : 

A = \V US \ = sin6 c = 0.2196 ± 0.0023. (3) 

The A parameter depends on the matrix element \V c b\ = 0.0395 ± 0.0017 (obtained 
from semileptonic decays of B hadrons |10]]) and on A: 

A= = 0.819 ±0.035. (4) 

y 

In order to constrain the parameters p and 77 without considering CP violation in the 
neutral Kaon system, three experimental input are used: 



2.1 



B° d oscillat 10ns. 



The mass difference Am^ between the mass eigenstates in the B® — B% system has been 
measured with high precision [JIT], [IT], . In the Standard Model it can be related to the 
CKM parameters in the following way: 



[l-pf + rf 



Am d 
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where mt is the the top pole mass scaled according to [15|| and z = mf/m 
F(z) is given by: 



F(z) 
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(5) 

The function 

(6) 



The values of all the parameters are given in table [I]. 

In the p — 7] plane, the measurement of Amj corresponds to a circumference centred 
in (1,0). By constraining 77 to zero, an evaluation of p can be obtained. Unfortunately the 



term fBdy B Bd , given by lattice QCD calculations, is known with a 20% order uncertainty 
|l~3| , |i"4|j and therefore it gives the largest contribution to the error on the p determination. 



2.2 B° oscillations. 



The mass difference Am s between the mass eigenstates in the B° s — system is expected 
to be much larger than Amd and in the Standard Model is related to the CKM parameters: 



P? + v 2 
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(7) 
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5 3692 ± 0020 GeV 
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Vb 


0.55 ±0.01 


Am d 


0.471 ±0.016 ps" 1 


Am s 


> 12.4 ps-^SS^o C.L.) 


\v ub \/\v cb \ 


0.093 ±0.016 


fB d \jB Bd 


0.201 ± 0.042 GeV 




1.14 ±0.08 



new input 



> 12.3 C.L.)p 
0.090 ±0.012 



191, 23 



0.215ig;8tgGeV [0, 0.210+ ' 039 



-0.032 



GeV p] 



1.14; 



-0.07 
-0.06 



Table 1: Physical parameters used in the formulae (5), (6), (7) and (8). In the second 
column the values used in [0] and in the third column the most recent values are given. 



Since the ratio 

^ _ Ib s V B Bs 
fB d \jB Bd 

is computed by lattice QCD with a better precision than the single terms, a measurement 
of could provide a much stronger constraint on the p — r] plane. However, given 
the very high frequency in the B° s — system oscillation, only a lower limit on Am s is 
available, as shown in £a&/e[I], which corresponds to a circular bound in the two parameter 
space or, if the assumption 77 = is made, to a lower limit for p. 



2.3 \V u b\ measurements from semileptonic b decay. 



Charmless semileptonic b decays have been used to measure \ V u b\ or the ratio \V u b\ / |K&|- 
The CLEO collaboration determined that parameter both by measuring the rate of lep- 
tons produced in B semileptonic decays beyond the charm end-point [I7| and from direct 
reconstruction of charmless B semileptonic decay [jTSfl . The two results are consistent but 
both methods are limited by theoretical uncertainties. In [ TDj] they are not combined and 
a value \ V u b\ / \V c b\ = 0.08 ±0.02 obtained from the former result is given. Recently CLEO 
gave a new result and an average with their previous results [W . At LEP, ALEPH [2^ 



L3 |2JJ and DELPHI [52] have measured the inclusive charmless semileptonic transitions 
b — > ulv. The average \V u b\ value from LEP measurements is given in |23| and the com- 
bination obtained using the latest CLEO result is given in table [I] in terms of \V u b\ / \ V c b\. 
The ratio of the two CKM matrix elements is related to the p and 77 parameters by: 



w, 



ub 



V, 



cb\ 



p2 _|_ ^2 



(9) 



and hence, in the p — r\ plane, the measurement of \V u b\ / \ V c b\ corresponds to a circum- 
ference centred at the origin. If 77 is assumed to be zero, it would be proportional to the 
p absolute value. 
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3 Data compatibility with a real CKM matrix hy- 
pothesis 

The assumption of a real CKM matrix implies that all the constraints described in the 
previous section are reduced to values of (or limits on) p. The compatibility of the obtained 
values can then be used to estimate the goodness of the assumption itself. In M and in the 
explicit reference to it in || , it is written that the hypothesis of a real CKM matrix can fit 
the data. However it is unclear which was the statistical approach followed in the second 
procedure in || to come to that statement [] and a completely different conclusion can 
be obtained with the same procedure and a standard statistical method. In addition the 
claim in the second procedure of ||, although using a different input-data-set, contradicts 
what is reported in 0. Assuming rj = , modifying accordingly eq. (5), (7) and (9) and 
using exactly the same input parameters as || (table |T] second column), the values 

p Am d = Q1 +0.18 (1Q) 

and 

p v ub = ± ( .42 ± 0.07) (11) 

are obtained from eq. (5) and (9), respectively. The limit on Am s has been obtained 
by means of the amplitude method |24j which allows to know the exclusion Confidence 
Level for any value of Am s . Therefore, by a convolution with the dominant uncertainty 
from the ratio £ in eq. (7) it is possible to obtain: 

p Am s > _ Q Q5 ( 12 ) 

at the 95% Confidence Level. 

In a naive approach on which the errors on p Arrid and p Vub are assumed to be uncorre- 
cted it is evident that the two values are fairly incompatible. The negative p Vub solution 
is clearly excluded by the p Ams limit and therefore it can be discarded. In order to in- 
clude all the correlations due to common terms contributing to the errors, namely \V c b\ 
and A, a two dimensional error matrix M has been written and the Best Linear Unbiased 
Estimator || has been used: 

„BLUE _ S f=l S j=lP»( M k? /,o\ 

with the variance ^ 

° 2p = E? =1 S? = i(M-i) y - (M) 
The error matrix M includes correlated and uncorrelated contributions: 



Me ■uncorr ^.uncorr , A A 

ij = d ij a i Uj +h a=1 A ai A aj (15) 

In the cited paper, a first procedure where both p and rj are fitted excludes the hypothesis rj = at 
the 99% C.L.. A second fit with ij forced to zero gives x 2 = 6.7 and the author concludes that this is 
compatible with a real CKM hypothesis but it is not specified why the mentioned % 2 value corresponds 
to a reasonable Confidence Level for such hypothesis. 
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where the indexes % and j run over the two p measurements and A a i is the change (with 
sign) on measurement % when the common systematic parameter a is moved by its error. 
The x 2 is obtained by: 

X 2 = S? =1 E| =1 ( 

With the ta&/e[l| (second column) parameters and taking the positive error in eq. (10), 
p BLUE = QM ± o 07 and x 2 = 4 3 ^ degree of f reec i om ) (17) 

are obtained. The values A a i and the more relevant uncorrelated contributions are given 
in table |2] and table ||| respectively. The correlation between the two measurements is 
small given the dominance of the uncertainty on fs d \jBs d in p Arrld as it can be deduced 
from the small off-diagonal term in the error matrix M (table f|). The corresponding x 2 
probability is 3.9% and this is clearly in contradiction with what stated after the second 
procedure of ||. 



parameter (a) 


A a pAm s 




v cb 

A 


4.3 
1.1 


-1.8 
-0.4 



Table 2: Variation (xlO 2 ) on p when the parameter a is moved by its positive error. 



Pi- P 



BLUE 



(m- 1 ; 



■ij 



pj - p 



BLUE 



(16) 



parameter 


i = p Amd 


i = p V u-b 


fB d \jB Bd 


+ 17.1 - 26.1 




m t 


±3.1 




Am d 


±1.7 




Vb 


±1.0 




\v ub \ 




±7.1 



Table 3: Uncorrelated contributions (xlO 2 ) to the p errors cij computed from the pa- 
rameters listed in the second column of table III. 





j = p Am d 


j = pVub 


% = p Amd 

i = p Vub 


3.3 x 10~ 2 
-0.8 x 10- 3 


-0.8 x 10- 3 
5.4 x 10~ 3 



Table 4: Elements of the error matrix M M{j computed from the parameters listed in 
the second column of table III. 

The effect of the limit on Am s is negligible even though the information contained 
in the amplitude is taken into account as suggested in ||: since the value of Am s that 
one would obtain by inserting p = p BLUE and 77 = in eq. (7) is Am s = 31 ps _1 , the 
present experimental sensitivity (13.8 ps _1 ) does not give any sizeable information on that 
p region. 

In order to take into account terms of the order up to 0(A 5 ) [25], the substitution 
p — > p = p(l — A 2 /2) must be done in eq. (5) with 77 = 0. 
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More recent values of fB d \ Bs d can be taken from WM (option a: fB d \ Bb, 



0.215±8:8m GeV) or from g§| (option b: f Bd ^B Bd = 0.210±g;^| GeV) where the lat- 
ter value has been obtained by considering the measured value of fo s and the theoretical 
evaluation of the ratio fB d /fD a - I n the following both cases will be considered. Including 
the last results on V uh [T^, ( see tableM, the incompatibility of the two p measurements 
is still present: 

-blue = Q 36 ± g Q5 and x 2 = 3 9 ^ degree of f reec i om ) (is) 

corresponding to a x 2 probability of 4.9% (option a) and 

-blue = Q 35 ± Q5 and x 2 = 4 3 ^ degree of free dom) (19) 

corresponding to a x 2 probability of 3.7% (option b). 

Since the dominant error in p Arrid is due to the lattice QCD computation, the hypoth- 



esis of a flat error with the same R.M.S. on fB d yBB d has been studied with a simple 
simulation. Several experimental results for p Amd and p Vub have been generated with cen- 



tral values equal to p BLUE . For p Arrid this is achieved by shifting the value of fB d yBB d 
and allowing it to vary within a flat distribution with the R.M.S. corresponding to the 
quoted error. All the other parameters of eq. (5) are allowed to vary with a gaussian 
distribution corresponding to their error. In each experiment the combined value and the 
X 2 are computed according to eq. (13) and (16), respectively. Looking at the x 2 distri- 
bution, it is possible to determine the fraction of simulated experiments with a x 2 higher 
than the value found in the evaluation with real data. This procedure has been repeated 
for the data sets of table [I] and for none of them the x 2 has been found to be higher than 
the experimental values of eq. (17), (18) and (19) in more than 5% of the cases. 

In it is mentioned that contributions from new physics can modify the value p Amd 
and the limit p Anis . In particular, being 5d, s the new physics contribution to Am^s, the 
changes on p Arrid <» depend on the fractional contributions Fd. s = 8d,s/ 'Am^s and with 
the data of the third column of table | Ap Amd ~ 0.5F d and Ap Am ° ~ 0.5(F d - F s ) are 
obtained. It is straightforward that in presence of a very large term Ap there cannot be 
any exclusion of the rj = hypothesis. It must be noticed, however, that those terms were 
not foreseen in the second procedure of || which has then to be compared with the present 
result from eq. (17). An evaluation of the allowed range of Fd from the data of ia&ZeH] 3 rd 
column option b), shows that small positive values of Fd are excluded (Fd > 0.04 at the 
95% Confidence Level) while all negative values are excluded except for the case of large 
positive F s (0(30%)) which would allow for negative values of p. This evaluation agrees 
with m. 



4 A method to determine p and rj 

To obtain p and fj without using the constraint on |efc| the eq. (5), (7) and (9) have to 
be used. It must be noticed that, in absence of a measurement of Am s , the number of 
constraints is such that a solution for the system of eq. (5) and eq. (9) is often found 
or, in other words, a 0(C) fit with p and fj as free parameters could converge to a pair 
of values such that the total x 2 function minimum is zero. The addition of terms taking 
into account the existing uncertitude in the value of the parameters entering into the two 



7 



expressions does not modify the number of degrees of freedom. With the present data 
(table |l] 3 rd column option b) such a solution exists and it is: 

p = 0.13ig;^ and fj = 0.38±g;^. (20) 

In the most recent search for oscillation |T2[ with the amplitude method an evidence 



for such a signal was not found but in the region 13 ps _1 < Am s < 17ps _1 the combined 
amplitude was greater than zero by more than 1.645 standard deviations. As explained 
J24|, the amplitude A and its error a a at any value of Am s can be related to the 
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log-likelihood referenced to its value for an infinite Am s (infinite oscillation frequency) 
by the expression A£°°(Am s ) = [1/2 — A(Am s )]/a\, Am n. From eq. (7) is possible to 
express Am s as function of p, fj, £ etc. and then the experimental amplitude and its 
error corresponding to a set of these parameter can be determined. Therefore a global 



log-likelihood function can be written as: £ = A£°°(p, 77, £) + £ Amd,Vub (p,ri, f Bd yB Bd , ..) 
where £^ m ^ v ^ j s the log-likelihood obtained from eq. (5) and (9) as function of the input 
parameters of those equations and their errors. Clearly £ !s - m ^ v ^ alone has a minimum in 
correspondence to the values of eq. (20) while the minimisation of £ gives: 

p = 0.14i{j;{jjj and fj = 0.37 ±0.05. (21) 

The addition of the term A£°° has the effect of excluding the p, fj space on the region 
corresponding to the excluded Am s values (i.e. negative p) but also of reducing consider- 
ably the errors of the two parameters on the opposite direction. The latter fact is related 
to the presence of a minimum on A£ oc (Am s ) for Am s = 14.75 ps _1 (see which is 
about 2.8 units below the asymptotic value for very high frequencies []26|| . However it has 
to be emphasised that, in absence of a clear signal for the B® oscillation, the presence 
of ghost minima on A£ QC (Am s ) cannot be excluded. Moreover, the minimum is slightly 
above the experimental sensitivity (14.3 ps _1 ). As a consequence, the results of eq. (21) 
must be taken as an indication of a method to include the Am s information (and of the 
relevant effects of its inclusion) rather than as a robust evaluation of p and fj with the 
present data. 



5 Conclusions 

The hypothesis of a real CKM matrix is tested on the basis of the present published 
and preliminary data and lattice QCD calculations. The Best Linear Unbiased Estimator 
has been used for the statistical approach. With all the input data used, included those 
suggested in ||, that hypothesis is excluded at more than 95% Confidence Level. This 
result agrees with the result indicated in with an older input data set and with a first 
procedure reported in || but contradicts statements in the same paper from a second 
procedure. A method to include the Am s information is proposed showing that a strong 
reduction on the error on p and fj is achievable. 
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